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Abstract 

We demonstrate the possibility of a self-consistent characterization of the photon-number statistics of 
a light field by using photoemissive detectors with internal gain simply endowed with linear input/output 
responses. The method can be applied to both microscopic and mesoscopic photon-number regimes. The 
detectors must operate in the linear range without need of photon-counting capabilities. 
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I. INTRODUCTION 



Although the measurement of photon-number statistics gives an incomplete description of the 
optical state, it can be helpful to provide fundamental information on the nature of any optical field 
and hence to discriminate between light of different kinds, either in the classical or non-classical 
domain. 

Measuring the photon-number statistics is a difficult task as it requires a complete knowledge of 
the entire detection process. In some situations, the evaluation of the first two statistical momenta 
can be enough for characterizing light statistics. In particular, mean value and variance, combined 
into the Fano factor, can discriminate between different statistics, even if sometimes the evaluation 
of higher order momenta is necessary. 

In this paper we demonstrate that it is possible to determine the photon-number statistics by 
means of a direct intensity measurement performed with linear photodetectors endowed with in- 
ternal gain. In order to bypass the problem of calibrating the detector, we devise a measurement 
technique that takes advantage of the linearity properties of the detector and is based on the mea- 
surement of the same state at different mean numbers of photons. In this way we simultaneously 
obtain the calibration of the detection chain and the reconstruction of the photoelectron statistics. 

Our technique has the additional advantage of being applicable to fields in the mesoscopic 
intensity regime (less than 1000 photons), which is scarcely explored. 

The competing techniques for determining the photon-number distributions are for instance: 

nnfi 

tomography with homodyne detection UJ, 121, 120, measurement with internal-gain photodetectors 
able to resolve peaks of different photoelectron numbers, plus analysis of peak integrals in order 
to reconstruct the photoelectron number distribution Drawbacks of the latter approach are: the 
low-limit in the mean number of detected photons due to the low quantum efficiency; the need of 
independent calibration of the detector; the need of a guess on the statistics to which the data must 
be fitted, as there is no direct indication from the measurement If the measurements are per- 
formed with linear photodiodes lacking internal gain, the drawbacks is the low-limit to the mean 
number of detected photons (p-i-n photodiodes can detect down to one thousand photons upon 
external amplification |5, ^). A further alternative is represented by ON/OFF measurements com- 
bined with maximum likelihood reconstruction [[TD. Section |II] presents the theoretical model we 
will use to interpret the experimental data, Sectionllllldescribes the experimental setup used for the 
measurements and Section |IV] presents the technique for analyzing the data and the experimental 
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FIG. 1: Sketch of the theoretical equivalent of a photodetector with rj < 1 endowed with internal gain, 
whose output signal is also externally processed to achieve an overall amplification by a factor a. 



results. Finally we draw some conclusions. 



II. THEORY 



We can model the detection process as a two step process: photodetection by the photocathode 
and amplification (both internal to the detector and external) (see Fig.[T]). 

We start from the observation that a real photodetector, having quantum efficiency i] < 1, 
behaves like a perfect detector (77 = 1) preceded by a beam spitter (BS) with transmittance T = 77 
[7, ^. The photodetection process can thus be described as the convolution of the photon statistics 
with the BernouUian action of the BS [8, ^. The statistics of the number of photoelectrons emitted 
by the photocathode, Pei{m), is thus linked to the statistics of the number of photons, Pph{n), by 



Peiim) = ^ B{m,n)Pphin) = ^ 



n 



m 



r^-{l-r^r^Pp^{n). 



(1) 



If we limit our analysis to the first two momenta of the distributions, the links between the statistics 
of photons and photoelectrons are given by 



m = ijn; a^i{m) = r] ap^{n) + ri{l — vi)n 



(2) 



where n = Yl'^=o Pph{n)n is the mean value and (y'^h{n) = Pphin){n — nY is the variance. 

The same notation is adopted for photoelectrons. 

A complete description of the amplification process requires the detailed knowledge, for each 
particular detector, of the amplification mechanism. However for the present discussion we limit 
ourselves to those detectors and to those operating regimes in which the amplification process can 
be simply considered as linear: we adopt a strongly simplifying approach in which the spread of 
the single photoelectron peak in the detector output is negligible as compared to its mean value. 
The relation linking the statistics of photoelectrons to that of the voltage outputs, even at the end 




FIG. 2: Sketch of the experimental setup. Left: Ught sources: (a) coherent state, (b) muhimode thermal state 
and (c) multimode pseudo-thermal state. D, rotating ground-glass diffuser; BBO, type I nonlinear crystal. 
Right: measuring apparatus: F, variable neutral density filter; MF, multimode fiber; PMT/HPD, detectors 
with internal gain; SGI, synchronous gated integrator; ADC, analog-to-digital converter. 



of the electronic chain external to the detector, is then 

Pout{v) = -Pei{am) , (3) 
a 

being a the amplification/conversion coefficient given by the overall amplification process and 
signal-processing electronics. Quite obviously, Eq. ^ shows that in general the distribution of the 
output voltage values is different from that of the photoelectrons. For the first two momenta of the 
distributions, the experimental final outputs are given by |19|] 

V = am] al^tiv) = a^al{m) , (4) 

where the symbols are defined as above. 

The most difficult task to be performed for reconstructing P^i from Pout is the experimental 
determination of the parameter a in Eq. This task can be tackled by following a procedure 
based on repeated measurements of the same state of light at different values of the quantum 
efficiency by inserting neutral filters in front of the detector. In particular we evaluate the following 
quantity 



Fy = — z — = = = arjP + ail — rj) = arjQ + a, (5) 

V arjn 

where F = a^f^ (n) /n is the Fano-factor that we have re- written as F = 1 + Q,Q being the Mandel 
Q-factor f.lQI . We observe that Fy is linearly dependent on rj and that a can be obtained as the 
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FIG. 3: Sketch of the experimental setup for the comparison between PMT and HPD. Nd:YLF, laser source; 
BS, beam splitter; F, variable neutral density filter; MF, multimode fiber; PMT/HPD, detectors with internal 
gain; SGI, synchronous gated integrator; ADC, analog-to-digital converter. 



limit of at ?7 — > 0. By multiplying and dividing by n we get 

Q 

= —V + a , (6) 

n 

which now expresses F^ as a function of the mean value of the measured output voltage. Note 
that V varies because of the variation of r] and not because of the change in the measured state 
(Q/n remains the same). The expression in Eq. ^ is general and all the information about the 
specific state under measurement is contained in the angular coefficient Q/n = {cr'^hin) — n)/n^, 
which is zero for Poissonian light, positive for classical super-Poissonian light and negative for 
nonclassical sub-Poissonian light. By measuring the same light state at different values of v and 
evaluating Fy, we can directly obtain the value of a and the ratio Q/n. Of course, different light 
states may have, for some values of the parameters, the same value of Q/n and this prevents the 
possibility of the unique determination of the light statistics. Nevertheless, if an assumption can 
be made on the statistics of the measured light, the method allows one to check it immediately. 

The expected expressions of F^ for the states presented in this paper will be given in SectionlTVl 
where the results are discussed. 



III. EXPERIMENTS 

In Fig. [2] we show the adopted experimental setup. The light from the specific source to be 
investigated is collected by a multimode fiber (100 /xm core diameter, OZ Optics, Canada). As the 
detector we used either a photomultiplier tube (PMT, 8850, Burle Industries, maximum quantum 
efficiency r] = 0.24 at 400 nm) or a hybrid photodiode module (HPD, H8236-40, Hamamatsu, 
Japan, maximum quantum efficiency r] = 0.40 at 550 nm). Both detectors are endowed with 
partial photon resolving capability and are linear over a wide range of intensities. The current 
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outputs of the detectors are integrated (SGI, Stanford Research Systems), sampled and digitized 
(ADC, National Instruments). The final outputs are then recorded by a computer. 

As a preliminary check of the performances of the detectors, we simultaneously measured the 
light at the two outputs of a 50% beam splitter (BS in Fig. [3]) dividing a ps-pulsed coherent beam at 
523 nm. The light was obtained from the second harmonics of a Nd:YLF continuous-wave mode- 
locked laser amplified at 500 Hz repetition rate (High-Q Laser Production). In Fig. |4] we show a 
typical pulse-height spectrum of the PMT (a) and of the HPD (b). The zero value of the output 
voltage V is set at the mean value of the response to dark of each detector, which was measured 
independently (see e.g. the horizontal scale in Figs.|4](a) and (b)). 

A standard procedure for analyzing a pulse-height- spectrum is to find a suitable fit obtained 
as a multiple convolution of the function best fitting the single-photon peak [4]. The procedure 
allows the reconstruction of a limited number of peaks (that depend on the resolving capabilities 
of the detectors) and fails in analyzing more intense fields when the pulse-height spectrum does 
not show well recognizable peaks. 

The pulse-height spectrum of the detectors lose the peak structure as soon as the field intensity 
becomes mesoscopic, still remaining proportional to the field intensity. 

In Fig. [5] we show several pulse height spectra recorded by the PMT (a) and by the HPD (b) 
in the same configuration of Fig. [3l in which neutral density filters are inserted on the beam path 
before the beam splitter. In Fig.[5](c) we plot the mean values of the output voltage as a function 
of the transmittance of the filters: the difference between the two curves reflects the difference in 
the overall quantum efficiency of the two detectors at 523 nm. The linearity will play a key role in 
the analysis procedure described in Section HVl 
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FIG. 4: Typical pulse-height spectrum measured with the PMT (a) and with the HPD (6) simultaneously 
recorded at the outputs of a 50% beam spUtter, dividing a ps-pulsed coherent beam at 523 nm. 



IV. RESULTS AND ANALYSIS 

We measured different light fields characterized by different statistics and at wavelengths that 
were chosen so as to match the highest quantum efficiencies of the detectors at best. The aim 
of the measurements is to identify the type of light field, to determine the overall amplifica- 
tion/conversion coefficient a and to reconstruct the statistics of the photoelectrons. The steps 
of the experimental procedure are: 1) to measure the same field at different attenuations of the 
neutral density filters (the actual position of the filters depends on the specific measurement); 2) to 
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FIG. 5: Pulse height spectra recorded by the PMT (a) and by the HPD (6) in the same configuration of Fig.|4] 

at different values of the transmittance of the neutral density filters, (c): Mean value of the output voltage as 
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a function of the filters' transmittance. The difference in the angular coefficient reflect the different overall 
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FIG. 6: Plot of the Fano factor, F^, for the output voltage as a function of the mean voltage v for different 
light states: pulsed coherent field at 523 nm (dots) and pulsed multimode thermal field at 420 nm (squares). 



evaluate the mean value and the variance of the output voltage; 3) to plot the experimental values 
of Fv (see Eq. ([5])) as a function of the mean value v; 4) to fit the plot to a straight line, thus 
obtain the value of the a factor, and to find, according to one of the theoretical models, the param- 
eters of the photoelectron distribution; 5) to use the fitting parameters to recover the photoelectron 
statistics. 

First of all we consider the measurements performed on ps-pulsed coherent light at 523 nm 
obtained as above (see Fig. [2] (a.))- For a coherent field we expect a Poissonian photon-number 
distribution 

Pph{n) = —r exp{-n) , (7) 
n! 

for which mean value and variance are n. Note also that the photoelectron distribution obtained 
by applying Eq. ([T]) to Eq. (|7]) remains Poissonian, so that Pei{m) = Pph{n) with the substitution 
n — i> m. The Fano factor for the final output voltage becomes 

F, = a, (8) 

independent of the mean value v. The experimental results obtained with the PMT are shown in 
Fig.[6](full circles). The linear fit gives a = (0.358 ± 0.002) V. 

The second kind of measurements is on the multimode thermal light given, as sketched in Fig.|2] 
(b), by a blue portion (420 nm) of the down conversion parametric fluorescence produced by a type- 
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I BBO crystal pumped by the third harmonics pulses (349 nm) of the same Nd:YLF laser (Si la]. 
The expected photon-number distribution is given by the convolution of /i independent thermal 
modes m 



+ 1)! 



Pphin) = -r-t ,.u-i , 1 J;../,-. , , (9) 



for which the mean value is n and the variance is a^^^ (n) = n {n/ jj, + 1). As in the case of coherent 
light, the photoelectron distribution obtained by applying Eq. ([T]) to Eq. Q remains multimode 
thermal and again Pei{m) = Pph{n). The Fano factor for the final output voltage turns out to be 

F, = - + a, (10) 

that is linear in the mean value. The experimental results obtained with the PMT are shown in 
Fig. [6] (squares). Their fit to Eq. ^ gives = 5.2 ± 0.1 and a = (0.356 ± 0.006) V. Note that, 
as these measurements were made in the same experimental conditions as those on the coherent 
field, the two values obtained for a are very similar to each other. 

We now use the fitting parameters to recover the photoelectron statistics. We divide the voltage 
output by the value of a to convert it to photoelectrons and then re-bin the obtained distributions 
in unitary bins. In Fig. |7] we show the photoelectron distributions, Pei.cxp for coherent light (a) 
and for multimode thermal light (6) reconstructed from some of the data sets used to obtain the 
calibration. The dots in the figure are the theoretical curves, Pgi (Poissonian and multimode ther- 
mal, respectively), evaluated at the measured mean values, by using the parameters, a or a and /i, 
obtained from the fits to Eqs. ^ and (flOl) , respectively. 

In the figures the values of the fidelity [11] of the reconstructed distributions. 



f=Yl \/ Pel,o.A^)Pelim) , (11) 

are also displayed indicating the good quality of the reconstruction. 

The same kind of measurements were repeated by using the HPD. To optimize the detection 
efficiency we exploited the light produced by the second harmonics of the Nd:YLF laser. Accord- 
ing to the setup depicted in Fig. [2](c), the laser light was passed through a rotating ground glass 
diffuser to obtain a pseudo-thermal speckle pattern [12]. The light, selected with a pin-hole so 
as to transmit a number of coherence areas, was then delivered to the detector by a multimode 
fiber. The photon-number statistics of the detected field follows the multimode thermal statis- 
tics in Eq. (HI). In Fig. [8] we plot the reconstructed photoelectron distribution for the multimode 
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FIG. 7: Reconstructed photoelectron distributions, Pe«,exp> for coherent light (a) and for multimode thermal 
light (6) reconstructed from some of the data sets used to obtain the calibration (bars). Symbols: theoretical 
curves, Pei, from Eqs. ([8]l and ( flOl ). 



pseudo-thermal field along with the theoretical curve evaluated from the fit of the Fano factor, 
a = (0.187 ± 0.002) V and /i = 3.9 ± 0.1 (see Inset of the Figure). Again the fidelity of the 
reconstruction is very good. 
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FIG. 8: Reconstructed photoelectron distribution, i-*ei,cxp> for pulsed multimode pseudo-thermal light at 
523 nm (bars). The displayed data set is one of those used to obtain the calibration. Symbols: theoretical 
curve, Pel, calculated from Eq. (ITOl ) with the number of modes (/i = 3.9) evaluated from the fit of the Fano 
factor (see Inset). Inset: plot of the Fano factor, F^, for the output voltage as a function of the mean voltage 



V. CONCLUSIONS 



We have demonstrated that it is possible to implement a self-consistent procedure to recover 
the distribution of photoelectrons that avoids calibration of the photodetectors. The procedure 
employs the evaluation of the Fano factor at different values of the mean photon numbers of the 
field to be characterized. The procedure has been satisfactorily tested on coherent and multimode 
thermal light. The next step will be to investigate less trivial optical states, such as those generated 
by the mixing of different fields jlS], in order to extend the validity of the method. 



12 



Aknowledgements 

The Authors thank Sergio Cova (PoUtecnico di Milano) for fruitful discussions. 
Present address of A. Agliati is: Quanta System, Solbiate Olona (VA), Italy. 



[1] M. Munroe, D. Boggavarapu, M. E. Anderson and M. G. Raymer, Phys. Rev. A 52 R924-R927 (1995). 
[2] Y. Zhang, K. Kasai and M. Watanabe, Opt. Lett. 27 1244-1246 (2002). 

[3] M. Raymer, M. Beck in Quantum states estimation, M.G.A. Paris and J. Rehacek Eds., Lect. Not. 

Phys. 649 (Springer, Berlin-Heidelberg, 2004). 
[4] G. Zambra, M. Bondani, A.S. Spinelli, F. Paleari and A. Andreoni, Rev. Sci. Instrum. 75 2762-2765 

(2004). 

[5] F. Paleari, A. Andreoni, G. Zambra and M. Bondani, Opt. Express 12, 2816-2824 (2004). 

[6] M. Bondani, A. AUevi, G. Zambra, M. G. A. Paris and A. Andreoni, Phys. Rev. A 76 013833 (2007). 

[7] G. Zambra, A. Andreoni, M. Bondani, M. Gramegna, M. Genovese, G. Brida, A. Rossi and M.G.A. 

Paris, Phys. Rev. Lett. 95 063602 (2005). 
[8] L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambrige University Press, New 

York, NY, 1995). 

[9] A. Aghati, M. Bondani, A. Andreoni, G. De Cillis and M.G.A. Paris, J. Opt. B 7 S652-S663 (2005). 
[10] L. Mandel, Opt. Lett. 4 205-207 (1979). 
[11] R. Jozsa, J. Mod. Opt. 41 2315-2323 (1994). 
[12] F. T. Arecchi, Phys. Rev. Lett. 15 912-916 (1965). 

[13] G. Zambra, A. AUevi, A. Andreoni and M. G. A. Paris, Int. J. Quantum Inf. 5 305-309. (2007). 



13 



